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Abstract 



(N 

■ The spin-1/2 Heisenberg model on an anisotropic triangular lattice 
is considered in the mean-field RVB approximation. The analytical es- 
timates for the critical temeperatures of the longitudinal s-RVB state 
(the upper one) and the 2D s-RVB state (the lower one) are obtained 

■ which fairly agree with the results of the previous numerical studies 
^> , on this system. Analytical formulae for the magnetic susceptibility 

and the magnetic contribution to the specific heat capacity in the lon- 
gitudinal s-RVB state are obtained which fairly reproduce the results 
of the numerical experiment concerning these physical quantities. 

1 Introduction 



The RVB state originally proposed by Pauling pQ for describing the structure 
of benzene molecule as a superposition (resonance) of its Kekule [2] , Dewar 
[3], and Claus [1] structures as shown in Fig. [TJhad been introduced into solid 
state physics by Anderson [5] in order to possibly explain the properties of 
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Figure 1: 7r-system of the benzene molecule. From left to right: two Kekule 
resonance structures, three Dewar resonance structures, and the Claus reso- 
nance structure. 



cuprate-based high-temperature superconductors. Since then this hypotheti- 
cal ground state is being sought in many materials. Those exhibiting frustra- 
tion of the exchange interactions of the magnetic momenta residing in them 
are widely suspected in forming a ground state of the RVB type. Among 
them the organic conductors of the family k-(BEDT-TTF) 2 Cu 2 XY [6] and 
the halocuprates of the formula Cs 2 CuCl 4 and Cs 2 CuBr 4 [7\ are considered 
as highly propbale candidates. Recently CuNCN phase had been obtained 
and a series of measurements had been performed of its spatial structure 
and magnetic susceptibility, electric resistivity, heat capacity (all vs. T) [8]. 
Although on the basis of analogy with other materials of the MNCN series 
(M = Mn, Fe, Co, Ni) QUI EEEJ one could expect more or less standard 
antiferromagnetic behavior, it turned out that in the low temperature phase 
the material at hand does not manifest any magnetic neutron scattering. 
The absence of the long range magnetic order (LRMO) (evanscence of the 
spin-spin correlation function) can be explained by the RVB character of 
the ground state of the system of the Cu 2+ local spins 1/2 in this material. 
Close inspection of the materials structure (Fig. 2 from Ref. [8]) reveals 
that each Cu 2+ ion can be effectively antiferromagnetically coupled to two 
its neighbors forming a chain while somewhat weaker antiferromagnetic cou- 
pling with four more neighbours from two adjacent parallel chains results in a 
Heisenberg model on an anisotropic triangular lattice with the Hamiltonian: 

^]Tj T S r S r+T (1) 

r r 

where the coupling vectors r take three values r^; i = 1 -j- 3; T\ — (1, 0); r 2 = 
(i,^);r 3 = (i,— with the interaction of the strength J along the lattice 
vector 7~i (two neighbors) and with a somewhat smaller strength J' along the 
lattice vectors r 2 and T3 (two neighbors along each). This is precisely the 
setting for which Hayashi and Ogata proposed that two (different) spin- 
singlet RVB (s-RVB) states are formed at different temperatures depending 
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Figure 2: ab planes of the CuNCN crystal. A stronger J\ extends in the 
a direction; somewhat weaker J 2 extends along the b ± a directions. The 
weakest J3 extends in the c direction and is not considered in the present 
paper. 



3 



on the amount of anisotropy y resulting in a temperature dependence of the 
magnetic susceptibility with characteristic discontinuities corresponding to 
installment of these two respective states [12J. 

In the present paper we reasses their numerical results analytically par- 
tially with use of the Ginzburg-Landau phase transition theory in order to 
further apply the obtained analytical formulae to the experimental data ob- 
tained for CuNCN [8]. 

2 RVB mean-field analysis of anisotropic tri- 
angular lattice system 

2.1 Equations of motion and self consistency equations 

Hayashi and Ogata [12] base their analysis of the Hamiltonian eq. ([I]) on 
returning to the electron representation from the spin representation by the 
standard formulae: ^ 

s i = ift^^Cifi, (2) 

where c^(cio-) are the electron creation (annihilation) operators subject to the 
Fermi anticommutation relations; cr Q/9 are the elements of the Pauli matrices 
and the summation over repeating indices is assumed. For the latter one can 
derive equations of motion based on the Heisenberg representation in which 
each operator obeys the following equation of motion: 

ihA = [A, H] (3) 

where [, ] stands for the commutator of the operators and '"' for the time 
derivative. Applying this to the creation and annihilation operators c^ a (c ra ) 
and performing commutation, mean field decoupling and Fourier transfor- 
mation as done in Appendix |A] results in mean field equations of motion for 
these operators : 

cos(kr)c+ k _ (T 

cos(kr)c_ k _ CT (4) 

These latter reduces to the set of 2 x 2 eigenvalue problems for each wave 
vector k: 



ihc^ = -- ^2 j -t£,t cos(kr)c k(7 -§ J2 T J -t^t 

T 
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( 6< A k \ / u k \ _ / u k \ 

Ui J k U J 

with 

£k = -3 2J Jrir cos(kr) 

r 

A k = 3 J r At cos(kr) (5) 

T 

(summation over r extends to ±7$; i = l-i-3) which results in the eigenvalues 
(excitation spectrum) of the form: 

£k = \J&+\^\ 2 (6) 

whose eigenvectors are combinations of the destruction and creation oper- 
ators with the above Bogoliubov transformation coefficients « k > v k . These 
equations result in the selfconsistency equations of the form: 

k / 

for six order parameters £ T , A T . 
2.2 Free energy 

Following Ref. [H] one can write immeditaly the free energy in terms of the 
above order parameters: 

(summation over r extends to ir^jz = 1t3). Minima of this expression 
correspond to various possible states of the system. This is more informative 
than studying the self consistency equations eq. (EJ) which may have more 
solutions. E.g. the trivial solution 

£t = A T = 

always exist although may be not of the lowest energy (see below). 
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Simplified model of s-RVB on the anisotropic 
triangular lattice 



The numerical analysis of Ref. [12] shows that in agreement with general 
theorems Ref. [13] the order parameters satisfy additional phase relations 
(see Appendix IA.3I for details) which allows to reduce the number of order 
parameters to only two: 

\/2Z = Z Tl ;V27i= |A T2 | = |A T3 | 

the first responsible for establishing the s-RVB state within the chains (longi- 
tudinal s-RVB) and the second for the same in the transversal direction (2D 
s-RVB). For the triangular anisotropic lattice as represented in the above 
approximation one can easily write the explicit expression for the free energy 
as relying on the general expression eq. It reads as follows: 

F = 6Je + 12JV ~ ^ / In (2 cosh ) d\ (9) 

BZ 

where 6 = ksT and 



El = 18 



J £ cos (k x ) + J rj I cos (— + -s-— ) + cos (- ^— ) (p 



It is easy to check that for — >■ the integral produces an expression 
which is uniform with respect to £ and rj so that the minimum of the free 
energy corresponds to both order parameters nonvanishing (see also below): 

£(0->O)^O;»j(0->O)^O. 

At high temperature we can use a high-temperature expansion: 

Ek\\ , n 1 / E k \ 2 l/£' k x! 



In 2 cosh [-=■ « In 2 + . 

V V 2 #// 2\2d J \2\26 

which allows to be integrated explicitly. This results in the Ginzburg-Landau 
approximate free energy -Fgl(£ ; Vi 6) which is of the fourth power in £ and rj. 
The critical temperatures in this approximation correspond to those where 
the free energy Fql of a state with one or two nonvanishing order parameters 
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becomes lower than that of some other state (e.g. of the paramagnetic state 
with vanishing order parameters). It happens at: 



Oc = I J, (11) 



where the longitudinal s-RVB state with: 



2y/26Jl-0/e c , s 

g = V 3J — —;v = o (12) 

installs (in perfect correspondence to the analytical result announced in Ref. 
[12]) and at 

°' c = 8(3 J - 2 J') (13) 

with 

«! = ^-2>l. 

e> c j> 

Below the lower critical temperature eq. f[T3"j) the 2D-RVB state installs with 
the order parameters: 




V = \/ht^ ( 14 ) 



which is in fair agreement with the numerical result of Ref. [12] in that sence 
that it shows some depletion of the longitudinal order parameter (£) below 
the critical temperature of the transition to the 2D-RVB state. The order pa- 
rameters' estimates are obtained as solutions for the minimum conditions for 
the Landau- Ginzburg free energy dF GL (^, rj, 0)/d(, = <9F GL (£, rj, 9)/di] = as 
functions of 9. The lower critical temperature eq. flT3|) is, as one can expect, 
an overestimate of the exact critical temperature usual for the mean field 
theory as combined with the high-temperature expansion. It is, neverthe- 
less, in a reasonable agreement with the numerical results of Ref. [12] where 
9'jOc was estimated to be 0.213 for y = 0.6 whereas our estimate is 1/3 for 
this amount of anisotropy. 
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For the longitudinal s-RVB phase eq. (IT2|) further analysis is possible. 
The integral term in the limit 9 — > reduces to 



TT I * 7T 

2 

with use of the integration method as described in Appendix O As combined 
with the term squared in £ in the free energy eq. (Q this yields the amplitude 
of the order parameter £ reached at the zero temperature: 

£o = 4-, (15) 

V Z7T 

which is in perfect agreement with the numerical result of Ref. [12J. 



4 Physical properties of the anisotropic tri- 
angular Heisenberg model 

Qualitatively the behavior of the anisotropic triangular Heisenberg model 
looks as follows: order parameters related to the intrachain and interchain 
interactions are both vanishing above the higher critical temperature 8 C . The 
excitation spectrum is gapless and there is no dispersion in either direction 
so that paramagnetic behavior is to be expected (see below). Below the first 
critical temperature (since — < 1 ) r\ = (longitudinal s-RVB state) the 
excitations have the energy spectrum: 

E k =3V2 |J£cos(k x )| (16) 

with no dispersion in the transversal direction although with the temperature 
dependent bandwidth. The excitations propagating with the wave vectors 
with the longitudinal components equal to ±| are gapless. Below the lower 
critical temperature one can expect that provided the ratio — is not too large 
the spectrum acquires a finite gap for whatever wave vector. Due to qual- 
itative changes in the spectrum of the system occurring at the two critical 
temperatures one can expect that the physical quantities which are predom- 
inantly contributed by the excitations with low energies will be significanly 
different in different phases. In what follows we shall derive analytical for- 
mulae for the magnetic contribution to the specific heat capacity and the 
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magnetic susceptibility in the longitudinal s-RVB state (between two critical 
temperatures eqs. fill I) . (1131) ) where only the order parameter £ is nonvanish- 
ing. In the temperature range of interest the properties are predetermined 
by the vicinity of the gap less excitations along the lines k x = ±| in the Bril- 
louin zone. However, the detailed form the spectrum must be immaterial for 
obtaining some semiquantitative estimates provided the qualitative features 
of the spectrum are captured. This can be done with use of a saw-like model 
dispersion law: 



3^2 



71 



;i7) 



which has the same gap behavior and the same bandwidth as the exact 
problem. It corresponds to the constant density of quasiparticle states in 
variance with the exact spectrum which is characterized by accumulation of 
the density at the top of the quasiparticle band. 

4.1 Specific heat capacity 

The magnetic contribution to the specific heat derives using the standard 
definition: 



where 







2kBT df ^ [/(Ek) ln /(Ek) + (1 " f{Ek)) ln(1 ~ f{Ek))] 



f(E) 



, E . ' 
exp(-) + 1 



is the Fermi distribution function. The specific heat capacity apparenly van- 
ishes in the gapless state (above C ). Replacing the summation by integration 
over BZ as described in AppendixOit becomes one per individual spin. After 
taking derivative with respect to T it becomes: 



V3 f E k (E k -TE'^')sech 2 (§) 



16tt 2 



BZ 



& 2 



d 2 k. 



where E^_ stands for the derivative of E^ with respect to £ and £' stands for 
derivative of £ with respect to T. 
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4.1.1 Saw-like spectrum 

For the saw-like spectrum in the area where the integration is to be performed 
the following 



9k 

E k (£ k - TE\£) = 18J 2 (1 - (£ " TO 

7T 



holds, so that 



where y = 2y / 3</£. The integration can be performed explicitly, which yields: 
/tanh(^) 14/ y \ tt 2 6 A9Li 2 (-e-e)\ 



V y \ ) 3y 2 

(18) 

For the divergent | one can expand hyprbolic tangent, the logarithm, and 
polylogarithm, and make sure that the divergent terms cancel each other 
and that the dominant contribution to the specific heat capacity is linear in 
temperature: 

C m = ^(l-T(\nyy)k B = -^k B 
3y 9V2J£ 

(correction to it is at best proportional to |e~^ which clearly disappears for 
the divergent |). On the other end of the interval - at the critical temperature 
8 C , by contract, ^ diverges according to the divergence law of £. Expanding 
the first bracket in the above expression eq. (|T8|) at small values of £ we 
obtain: 

C m = k B -^(l-T(\nO'). (19) 

This exact result can be combined with the classical expression for the tem- 
perature dependence of the order parameter in the vicinity of the critical 
point (for the reasons which will be clear below we use here a generalized form 
of the classical temperature dependency of the order parameter (1 — QjQ^I 2 
which shows the same critical exponent as the classical one): 



f = foy'l - (W 

n = W 7T=kf (20) 
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First term in the brackets eq. f ll9p when multiplied by £ 2 is small in a vicinity 
of 9 C . The second term yields: 

= a fey i a m a /CoV 

e 2 U; 2 U; U/ 

so that, combining this with eq. (Tl9|) we obtain the finite jump of the mag- 
netic contribution to the specific heat capacity in the critical point: 

a ^ , n \ , 3J 2 a ~ 16a, . . 

AC m (0c) = k B -p-tl = ^k B . (21) 

This finite jump is qualitatively in agreement with general theory Precise 
agreement with the numerical resutls of p2] is not achieved due to approx- 
imate form both of the spectrum and the temperature dependence of the 
order parameter used above. On the other hand required agreement can be 
achieved by fitting the value of a (see below). 



4.1.2 Exact spectrum 

If the exact spectrum is inserted in the general expression for the magnetic 
contribution to the specific heat capacity we get: 

C m = ks^— 1 — — - f cos 2 x sech 2 (^—^ — 1 dx, 



tt0 2 Jo V 26 

We could not find analytical expression for this integral. However, the mag- 
nitude of the jump of the magnetic contribution to the specific heat capacity 
in the critical point can be estimated from an expansion of the above integral 
for small y (£). The result in the vicinity of the critical point reads: 

Combining this with the above expression for the logarithmic derivative of 
the order parameter we obtain: 

8a 

AC m (O c ) = —k B - (22) 

Analysis of the low temperature limit for the exact heat capacity will be 
performed elsewhere. 
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4.2 Magnetic susceptibility 

We use the standard definition of the magnetic susceptibily per spin: 

k K 

The summation replaces by integration which yields: 

BZ 

Following the recipe of Appendix IB. II we get for the susceptibility: 



4n 2 6 

BZ 



which is analyzed below. 



4.2.1 High temperature (9 > 9 C ) gapless state 

In this state the argument of the function under integral equals to zero for all 
values of k. Thus inserting = to the derivative of the Fermi distirbution 
function 

df(E) = exp(f) 
dE ^exp(f) + l] 2 

we obtain: 

in perfect agreement with the numerical result reported by Hayashi and 
Ogata [12]. 



4.2.2 Longitudinal s-RVB state 

In this state the gap is absent along the lines k x = ±| in the reciprocal 
space. Their vicinity provides major contribution to the susceptibility. It 
will be estimanted below for the saw-like approximation of the spectrum and 
for the exact one. 
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Saw like approximation of the spectrum Performing the integration 
as described in Appendix O we get for the susceptibility: 



n=l 



Summation of the constant term yields | according to Appendix IB . 1 1 and the 
rest is just the geometric series so that the susceptibility becomes: 



X = ~F7 



/4V2 [ 1 ex P (- 



3V2JA 
9 J 



which transforms into 



J ^ ! I 2 l + ffl q,(-2jf2) 



y-^^tanh ^ 



Qualitative behavior of the saw-spectrum model susceptibility The 

temperature behavior of the susceptibility in the saw approximation depends 
on that of £(#). If the latter flows to a finite value while 6 — >■ the susceptibily 
remains finite as well in agreement with finding of Ref. [T2] for the longitudi- 
nal state. The limiting value equals which is significantly higher than the 
numerical limiting value obtained in Ref. (T2]. By contrast, if — > a > 
(the case of the temperature dependence provided by the high-temperature 
expansion eq. f|T2|) ). a paramagnetic behavior can be expected, although 
with the number of momenta effectively suppressed by hyperbolic tangens 
which is always less than unity. 



Exact susceptibility in the longitudinal s-RVB state Using the inte- 
gration trick as described in Appendix [C] and inserting the exact spectrum 
of the longitudinal state and noticing the argument of the absolute value is 
positive in the integration limits we obtain for the integral: 

* exp (^) dkx = I exp I ^ - I dfc. 



2 




L 



n 



3V2JM 
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where Jo is the modified Bessel function and L is the Struve function of the 
0-th index. So that we arrive to the expression for the susceptibility in the 
longitudinal s-RVB state 



m = Mq,)«„ ^ pb*«) _ l0 p^i)) , (25) 

which is exact. 

Qualitative behavior of the exact susceptibility The differences of 
the modified Bessel function and the Struve function in eq. (1251) divided by 
9 have finite limit while 9 — > which equals to 

37rJ£n 

Inserting it into the series eq. ( 12 5 p and performing summation according 
to Appendix |B] we obtain that at zero temperature the susceptibility in the 
longitudinal s-RVB state has finite limit 



X 



/4\/2 2/4 



37tJ£ 3J 

This is in perfect agreement with the numerical result of Ref. [12]. The ratio 

^exact^ = o) _ 2 i 

which is expectable due to accumulation of the quasiparticle states at the top 
of the exact spectrum vs uniform density of states for the saw-like spectrum. 
On the other hand in the vicinity of the critial temperature 9 C where order 
parameter is small in eq. ( |25|) the series terms have finite limit of unity 
independently of the n value, which as combined with Appendix IB. II yields: 



so that 

X (9 = 0) 2V29 C 



X (9 C ) 3ttJ£ 



< 1. 



o 



The susceptibility at the critical temperature 9 C is continous as well in agree- 
ment with results of Ref. [12] as they are with the results of the saw- 
approximation for the spectrum. 
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High-temperature estimate for the temperature dependence of the 
susceptibility in the longitudinal s-RVB state At high temperatures 
(in the vicinity of the upper critical temperature) the temperature depen- 
dent bandwidths of the quasiparticle spectrum scaling the arguments of the 
functions in the series eq. ( 125]) are small. Nevertheless the expansions of the 
modified Bessel function and the Struve function at small values of argument 
cannot be applied since the summation extends to arbitrary large values of n. 
For that reason in order to derive temperature dependence of the susceptibil- 
ity in the vicinity of the upper critical point we apply the high-temperature 
expansion for the energy derivative of the Fermi distribution eq. (I24|) as 
described in Appendix IB.2I and arrive to the formula: 



/4 

29 



4 E 



n=l 



y_ 

20 



(2 



n+2 



1)B. 



n+2 



(fO (n + 2) 



(26) 



to be used for numerical estimates. 



Asymptotic estimate for the temperature dependence of the sus- 
ceptibility in the longitudinal s-RVB state In order to eluciate the 
details of the temeperature behavior of the susceptibility in the longitudinal 
s-RVB state close to the zero temperature (provided this state exist in this 
temperature region) and eventually valid above the lower critical value we 
peformed analysis using asympthotic expansions for the Bessel and Struve 
functions at large values of their arguments (low temperatures) as described 
in Appendix IB. 31 and finally arrived to: 

^^(H'jlK-fp-,)^ (27) 

The performed analysis is valid provided the ratio ^p- has nonvanishing limit 
or diverges while 9 — > 0. In this case the arguments of the Bessel function in 
eq. (|25p remain large and the asymptotic expansion used in Appendix IB. 31 
on which the final result is based remains valid. 

4.2.3 Two-dimensional s-RVB state 

In the region below the lower critical temperature 9' c both order parameters 
£ and rj of the simplified model are nonvanishing and the spectrum has a 
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finite gap for all values of the wave vector in the BZ. Thus the susceptibility 
can be fairly approximated by a simple expression: 

_ 2/i| exp(-f ) 
X [exp(-f) + l] 2 

where A is some characteristic value of average gap (constant in the BZ) effec- 
tively taking into account both the gap and the dispersion of the quasiparti- 
cles. This expression describes exponentially decaying susceptibility provided 
— > oo while 9 — > 0. This is in qualitative agreement with the numerical 
results of Ref. [12]. Further improvements can be expected from taking into 
account temperature dependence of the average gap below the lower critical 
temperature due to expected temperature dependence of the transversal and 
longitudinal order parameters. 



5 Comparison with numerical experiment 

It has been mentioned above that the precise temperature behavior of the 
physical quantities as coming from the analysis of the exact spectrum or 
with a saw-like model depend on the details of the temperature dependence 
of the order parameters £ and 77. The numerical results of Ref. [12] show 
that the both order parameters remain finite at 9 = 0. The high-temperature 
expansion as yielding the critical temperatures for the emergence of the both 
order parameters yields a linear evanescence of £ in this limit in contradiction 
with results of numerical experiment and with the exact result eq. (I15p . 
Obviously the high-temparature expansion fails at low temperatures. On the 
other hand close inspection of the numerical results of Ref. [121 US] shows 
that in addition to the vertical tangent characteristic to the £ (9) dependency 
at the upper critical temperature one has to count to visibly horisonthal 
tangent at zero temperature. This indicates that the exponent a in eq. f l2"Uj) 
must be larger than unity. With these precautions we shall establish relation 
between our analytical results and those of the numerical modeling in Ref. 

5.1 Heat capacity 

According to calculations of Ref. [12] the magnetic contribution to the spe- 
cific heat capacity experience a jump of 2fcg in the upper critical point. 
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C m /k B 




0.05 0.10 0.15 0.20 0.25 0.30 0.35 



Figure 3: C m vs 8 as obtained by numerical integration with the exact spec- 
trum and with generalized classical approximation for (a = vr 2 /4). 

Combining this with our result eq. (122]) we can derive a = tt 2 /4 « 2.46. 
Corresponding graph represented in Fig. |3] shows perfect agreement with the 
numerical results of Ref . p2] . 

5.2 Magnetic susceptibility 

The results of the model calculation on magnetic susceptibility per spin are 
represented in Fig. H]for the generalized classical approximation for £(8). 

The temperature dependence of the suceptibility outside the longitudinal 
RVB region (8 < 8' c , 8 > 8 C ) is given respectively by the average gap model 
and by paramegnetic susceptibility. Within the longitudinal RVB region 
{8 > 8' c , 8 < 8 C ) the asymptotic expansion eq. (12?]) is used in the vicinity of 
the lower critical temperature and the high-temperature expansion eq. ( 1261) 
in the vicinity of the upper critical point. The picture fairly shows the key 
features of the temperature dependence of the magnetic susceptibility of the 
longitudinal s-RVB state as they appear from the numerical calculation of 
Ref. [12]. The low-temperature limit is almost achieved pretty above the fi- 
nal exponential drop to zero characteristic for the low temperature 2D-RVB 
phase. The asymptotitc and high-temperature expansions for the susceptibil- 
ity of the longitudinal s-RVB state derived for low temperatures remain fairly 



17 



xly-B 2 
1.4 - 




0.0 0.1 0.2 0.3 0.4 0.5 

Figure 4: \ vs 9 with generalized classical approximation for £(T) (a = 
7r 2 /4). and in . The susceptibility in the low-temperature phase (2D-RVB 
with both longitudinal and transversal order parameters nonvanishing) is 
drawn in the average gap approximation with the exact magnitude of the 
latter adjusted to assure continuous behavior of the susceptibility at the 
lower critical temperature. The interval of 9/ J between 9' c and 2- J is drawn 
according to eq. fl27|) with k max = 1 . The interval of 9/ J between ^j- J and 9 C 
is drawn according to eq. fl26|) up to 14-th power in reciprocal temperature. 



18 



valid (on respective sides) almost until quite arbitrary temperature ^ J with 
the acceptable amount of discontinuity between two different expansions. 

6 Conclusion 

In the present paper we succeeded in obtaining analytical estimates for the 
critical temperatures of the installment of the longitudinal and 2D s-RVB 
states in the spin-1/2 Heisenberg model on the anisotropic triangular lattice 
in the Ginzburg-Landau approximation for arbitrary value of the anisotropy 
(in the assumption of the existence of the both longitudinal and 2D states). 
These estimates are in a fair agreement with the numerical results of Ref. 
[T2] related to the same system, but apply for the arbitrary anisotropy value 
y in the allowed range. The method does not allow for obtaining correct 
estimates of the temperature dependence of the equilibrium values of the 
order parameters of the system. That of the longitudinal s-RVB state was 
interpolated with use of the generalized classical temperature dependence 
corresponding to the critical exponent 1/2 and scaled so that its value at zero 
temperature equals to that obtained by minimization of the zero temperature 
free energy. 

For the longitudinal s-RVB state of the spin-1/2 Heisenberg model on 
the anisotropic triangular lattice, existing between the upper and lower criti- 
cal temperatures the analytical estimate of the magnetic contribution to the 
specific heat capacity is obtained with use of the saw-like approximation for 
the quasiparticle spectrum. The low-temperature behavior of the magnetic 
contribution to the specific heat capacity coincides with that obtained nu- 
merically in Ref. [12]: it has a finite slope at 6 — > 0. On the other hand the 
magnitude of the abrupt jump of the magnetic specific heat capacity in the 
paramagnetic state (above the upper critical temperature) had been derived 
for the exact spectrum in agreement with general theory. In the generalized 
classical approximation of the temperature dependence of the order param- 
eter its magnitude (2/cb) known from numerical experiment Ref. [12] can be 
reproduced under setting the a parameter of the generalized classical formula 
equal to 7r 2 /4. 

For the longitudinal s-RVB state of the spin-1/2 Heisenberg model on 
the anisotropic triangular lattice analytical expressions are obtained for the 
magnetic suceptibility for the saw-like model of the quasiparticle spectrum as 
well as for the exact spectrum. (In the latter case - the asymptotic estimates 
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have been obtained valid respectively in the vicinity of the lower and higher 
critical temperatures). Under assumption of the generalized classical tem- 
perature dependence of the longitudinal order parameter as the most realistic 
from the point of view of reproducing the results of numerical experiment 
on the magnetic contribution to the specific heat capacity the temperature 
dependence of the magnetic susceptibility in the longitudinal s-RVB state 
are fairly reproduced: the finite limit at the zero temperature (under con- 
dition of the existence of the longitudinal s-RVB state at this temperature) 
is obtained as well as the decrease of the susceptibitity in the longitudinal 
phase as compared to the maximal value achieved at the upper critical tem- 
perature, separating the longitudinal s-RVB phase from the paramegnetic 
one. Narrow problematic region where both the low-temperature asymptotic 
and the high-temperature expansions are not good enough is found around 
9 « 3J/11. 
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A Derivation of equations of motions 

A.l Equations of motion for Fermi operators 

The equation of motion for the Fermi operators as coming from the Heisen- 
berg spin Hamiltonian eq. ([[]) are based on the commutation relations be- 
tween the Fermi operators and the spin operators as given by eq. fl2]). In- 
serting the explicit form of the tensor components of the spin operator: 

S i = c ia Cip; = cfpCi a ; 

S i = (<&Cia- 4 p Cip), '2 

one obtains: 

SjSj] — {cipSj + c ia Sj)/2 

[Ci^SjSj] = (Sfc ia - C i/3 Sj)/2 
with this we obtain in terms of the Fermi operators: 

[ c iai SjSj] = (CipCjpCj a Cict(Cj a Cj a C^gCj / g)/2)/2 

[c%p> SjSj] = (cj Q Cjpc ia — Cip(cj a Cj a — CjpCjp)/2)/2 
which must be inserted for each terms in eq. (Q. 

A. 2 Mean field decoupling of equations of motion 

The terms contributing to the right hand side of equations of motion for 
Fermi operators as derived in the previous Section are all of the third overall 
power with respect to these operators. If one is interested in having an option 
for superconducting (anomalous) state the mean-field decoupling has to read 
as: 

cfc 2 c 3 =3> (cfc 2 ) c 3 - (cfc 3 ) c 2 + cf (c 2 c 3 ) 

where the superconducting anomalous averages appear. In the RVB state we 
are going to take care of it is reasonable to set 

In order to assure the ^conservation we set 

(Ci a Cj a ) = (CipCjp) = 
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<4c- CT ) = 

After further assuming the translational invariance and effectively introduc- 
ing order parameters £ T : 

( C ta C ja) = { C % C jP) = 
£t = ( C r+To- C rcr) 



(c r +TQ C r/?) — (c r aC r + Tj g) — A T — A_ T 
( c r+r/3 c r«) = ( c r^ c r+r«) = ~A T = ~^-t 

Inserting these results in the equation of motion and using the Fourier trans- 
forms for the annihilation operators 



c k(T = — L ^ exp(-zkr)c r 
v A 



(A is the number of sites in the crystal) and hermitean conjugate one for the 

creation operators we arrive to: 

the final equations of motion of the form: 

ihc ka = -- ^2 J t£t cos(kr)c kQ +§ ]T r J T A T cos(kr)c+ k/3 



ihckp = -r Jt£t cos(kr)c k/3 -§ ^ T J T A T cos(kr) 



where summation over r is extended to ±7*; i — 1 -r- 3 



C -ka 



A. 3 Phase relations for order parameters 



As it can be seen from the definitions eq. (jSJ) of the dispersion functions 
and the requirement that the final form of the expression eq. ([6]) must be 
the square root of a sum of square terms referring to each r as is mutatu 
mutandis eq. (12) of Ref . [T5] : 



El 



f'2,2 2,i r/2 2 f 2/^, kjy\/3 

J £ cos (k x .) + J 7? ( cos (— + ) 



i 2/^ 

+ cos ( y 



k„V^, 
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where we set £ = V2D 1 and rj = \^2D 2 3. It can be achieved if the cross 
terms between the contributions with from different r's cancel each other. 
This can be done by specific selection of the relative phases and amplitudes 
of the six order parameters £ T and A T . We assume that 

Then the condition for disapperance of the cross terms in the expression for 
E%. reads: 

£t£t' = -VtVt' cos(y? r - yv) 
for all unordered pairs r, r' (three conditions). They can be nontrivially 
satisfied e.g. by the following choice (not unique due to SU(2) symmetry of 
the problem): 

Vi = 0;6t^0 
V2 ^ 0;£ 2 = 
Vs + 0;£ 3 = 

7T 

- = y?2 - V^3 
& = 2£ 2 

2 2 o 2 

V2 = V3= 2 V 

which yields the required form of the excitation spectrum in the RVB phases. 

B Series summation 

In the present work we use asymptotic series as well as series not uniformly 
convergent. The sums of the latter are undestood as limits of the functions 
they represent on the border of the convergency range. Examples follow. 

B.l Using geometric series for summation of divergent 
series 

The Fermi distribution can be represented as a shifted geometric series: 

-ffH" - i5a=-H)£.-»--(-¥) 

= £(-1) exp^ _j 

n=0 V 7 
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which is the same as 

oo 

= x-x 2 + x 3 ... = Y(-l) n x n 
1 + x ^ 

n=0 

with 



x = exp 



E 
J 



This series converges for all \x\ < 1 or 6 > 0, but at x — 1 it becomes a 
divergent series 

oo 

^(-l) n = 1-1 + 1-1 + 1... 

n=0 

The function on the left hand, which defines the series, is continous in x — 1 
and the the value of this function is taken as that of the sum of the divergent 
series: 

E<-i)"4 

n=0 

Analogous consideration applies to the derivative of the Fermi distribution 
with respect to the argument of the exponent: 

exp(f) ~ / (n + l)£ 



l2 

ra=0 



£ (-!)•(„ +1) exp (-^^) 

■n—n \ / 



[exp(f) + l]' 

Shifting the summation index by unity makes it a series: 



x 



- ^(-l) n+1 (n + l)x n+1 = - ^(-l) n nx n , 



(1 + x) 2 

V / n=0 n=l 

which is the derivative of the geometric series: 

1 oo 

= l _ x + x 2 ... = V(-l)"x n 

1 + X „ 
1 V I A 



1 + X 



)/ ^ CXJ CXJ 

[i- + x) n „ 
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multiplied by x. The series on the right hand side apparently diverges at 
x = 1 where it becomes: 



- ^{-l) n n = -1 + 2 - 3 + 4... 

n=l 

However, its sum, being defined as a limiting value of the function on left 
hand side, must be: 



oo 1 

-£(-!)•» = ; 

which is used throughout. 
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B.2 High-temperature expansion for the susceptibility 
in the longintudinal s-RVB state 

In order to obtain the high-temperature expansion for the susceptibility we 
notice that the derivative of the Fermi distribution eq. ( |24|) upto a numerical 
coefficient coincides with the generating function for the Euler polynomials 
EJx) AS23.1.10: 



g/_ v /l\" ti K +1 (0) 
de ^\9J 2n\ 

n=0 v 7 



-e n , 



which is a good series expansion for small values of e/6 (higher tempera- 
tures/small quasiparticle bandwidth). Inserting this and the exact quasipar- 
ticle spectrum e(k x ) of the longitudinal s-RVB state in the definition of the 
magnetic susceptibility eq. f l23|) we obtain with account of the integration 
trick over the Brillouine zone in the longitudinal s-RVB state Appendix O 
for the susceptibility: 



a + i(q) /■* B . 

n=0 v 7 



Inserting the values of integrals which are ^^i 2 ^^ and employing the equal- 
ities for the factorials and Gamma functions AS6.1.18 we arrive to: 



x = /4 



n=0 



n+1 



y\n E n+1 (0) 



2J [(n/2)! 



2 ' 



hereinafter ASnn.n.nn stands for the formula nn.n.nn of Ref. 1171. 
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Then using AS23.1.20 we get 



A 29 



1 + 4 



£(!) 



y\n (2 



n+2 



n+2 



-\29J [(n/2)!] 2 (n + 2)_ 
(here £> n are the Bernoulli numbers) to be used for numerical estimates. 

B.3 Asymptotic low-temperature expansion of the sus- 
ceptibility in the longintudinal s-RVB state 

We make use of the asymptotic series AS12.2.6for the difference of the Bessel 
and Struve functions: 

I (x)-L (x)~-}^ — 



*t^m-k) (D 



2k+l 



valid for large values of argument: 



i. e. for the lower temperatures. It is remarkable that the expansion behaves 
better for larger values of n. Inserting this expansion in the series for the 
susceptibility we get: 



X 



7l9 



n=l k=0 



ni-k) (n) 



2fc+l ' 



which results in: 



2/i 



2fc+l 



The summation over n then performs with use of AS23.2.19 

oo 

J2(-lT +1 n- k ={l-2 1 - k )((k), 



n=l 



where ((h) is the Riemann zeta function, which results in 

X 



2/4 c 
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where the term with k = yields: 

E, i\n+l 2^ 
(-1) n— = -, 

n=l ^ 2/ 

according to Section [Bj] so that: 



X 



where the summation index k is shifted by unity. Next following AS23.2.16 

we get 

v 2& | /4 f - (-i) fc r(| + k)(2*y» (2e\ 2k+1 , k 

where B 2 k is the corresponding Bernoulli number. 

Further simplification is performed with use of the formulae for the Gamma 
function. Using AS6.1.17: 

r(z)r(i-z) 



sin(7rz) 



for z = \ + k we obtain 



2 ' ' r(| + k) sin(7r(i + k)) 



l) k T{\ + k) 2 



so that 

X = — + ( 7 j (1-2 

Inserting the definition of factorial in terms of Gamma function: (2k)\ 
T(2k + 1) and applying the duplication formula AS6.1.18 

(2k)\ = ^=2 2k+1 T{ l - + k)T{k + l), 
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we get: 

2,4/ i f- r(i + <Q (*e\ 2t 2k A 

to be used for numerical estimates. 

C Integration in the Brillouin zone of the 
trigonal lattice 

The Brillouin zone (BZ) for the triangular lattice is a symmetric hexagon 
with a side An/ 3 with the diagonal parallel to the k x axis in the reciprocal 
space. Its area equals to 8n 2 / \/3 as employed in the expressions throughout 
the paper. The integration over BZ falls into three parts: 

• one (1) over the bulk of the BZ: $\dk x Jj<L dk y 

2 v/3 

2tt 12L 

• one (2) over stripes 2 x f K 3 dk x Jj^ dk y 

2 \/3 

• one (3) over triangles 4 x dk x j^ dk y 

3 

We notice, however, that in the longitudinal s-RVB state the excitation spec- 
trum is an even function of k x which consists of two branches each symmetric 
with respect to reflection in the straight lines k^. = ±7r which correspond to 
the ridges of the respective branches of the excitation spectrum. For that 
reason while integrating any function independent on k y over k y over the 
flaps of the BZ (the values of ^ < |k x | < -y - triangmes integration) we 
see that for each given constant value of the integrand the length of the 
integration range corresponding to this value on one side of the ridge as 
summed with that on the other side of the ridge yields precisely ^= which 
is the width of the BZ in the k y direction. For that reason the integration 
over k y in the flaps of the BZ is equivalent to integration over the rectangle 
< kj, < 7r; 1^ I < which after completing it with the integration over 
the rectangle | < k x < 4p |k y | < (stripe integration) becomes equal to 
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the integration over the bulk of the BZ, Thus the integral over the entire BZ 
in the longitudinal RBV state satisfies the condition: 

I — — 

as used in our calculations. 
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